Domain structure of superconducting ferromagnets 
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In superconducting ferromagnets the equilibrium domain structure is absent in the Meissner state, 
but appears in the spontaneous vortex phase (the mixed state in zero external magnetic field), though 
with a period, which can essentially exceed that in normal ferromagnets. Metastable domain walls 
are possible even in the Meissner state. The domain walls create magnetostatic fields near the 
sample surface, which can be used for experimental detection of domain walls. 



PACS numbers: 74.25.Ha, 74.90. +n, 75.60.-d 

Recently there has been a growing interest to materi- 
als, in which superconductivity and ferromagnetism co- 
exist [0 |H). A number of unusual phenomena and struc- 
tures have been predicted and observed, spontaneous vor- 
tex phase as an example |0fl. But the theory mostly 
addressed macroscopically uniform structures, whereas 
ferromagnetic materials, even ideally uniform, inevitably 
have a domain structure, which is a ground-state prop- 
erty of ferromagnets. So a further progress in studying 
materials with coexisting ferromagnetism and supercon- 
ductivity requires an analysis of the domain structure. 
The present work is the first step in this direction. 

An object of the study is a material, in which the mag- 
netic transition occurs earlier, i.e. at a higher tempera- 
ture, than the superconductivity onset. This was called 
"superconducting ferromagnet" , in contrast to "ferro- 
magnetic (or magnetic) superconductors" where the su- 
perconductivity sets in before the magnetic transition, 
which have been studied mostly in the past || . Competi- 
tion of ferromagnetism and superconductivity may result 
in various structures with the magnetic moment rotating 
in space (spiral structures, cryptoferromagnetism and so 
on) . This also can be considered as a "domain structure" , 
but with a period determined by intrinsic properties of 
materials. However, our goal is the domain structure 
due to magnetostatic fields generated by nonzero aver- 
age bulk magnetization M. In this case the domain size 
depends on a sample size. We shall consider type-II su- 
perconductivity, bearing in mind ruthenocuprates JjiRE) , 
which are type II high-T c superconductors. 

Before analyzing the domain structure it is useful to 
summarize the magnetic properties of a single-domain 
superconducting ferromagnet. The total free energy of 
the superconducting ferromagnet can be written as H 



F(M, B) = f E + K 
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where A is the London penetration depth and B is the 
magnetic induction. The energy /e(M,VM) is the ex- 
change energy, which depends on the absolute value of 
M and on gradients of M. As a rule fl, in magnetic 
materials this is the largest energy, which fixes M. The 
anisotropy energy K(M /M) is smaller and depends on 



the direction of M. We shall consider a stripe magnetic 
structure, which is possible only if K essentially exceeds 
the magnetostatic energy ~ M 2 0. The latter is deter- 
mined by the magnetic field H — B — AttM [the third 
term in Eq. (jj])]. The expression Eq. (|l|) includes also 
the kinetic energy related to the superconducting current 
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where $o is the magnetic-flux quantum, ip is the phase 
of the superconducting order parameter, and the vector 
potential A determines the magnetic induction B = V x 
A. The kinetic energy of superconducting currents is 
absent in a normal ferromagnet. 

Minimization of the energy with respect to the vector 
potential A yields the Maxwell equation 
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V X (B - 4ttM) = VxB 



Together with the equation 

V x j s = - 
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this yields the London equation which determines B: 

A 2 V x [V x B] + B = . (5) 

Here we took into account that V x M = inside do- 
mains. In contrast to Ref. jj, we neglect the differential 
susceptibility (M does not depend on a magnetic field), 
which renormalizes the London penetration depth. 

These equations and the boundary conditions at the 
sample boundary (continuity of the tangential compo- 
nent of H and of the normal component of B) yield the 
distribution of B and H . This distribution is shown in 
Fig. for the case of M parallel to the sample bound- 
ary and for zero external magnetic field. The magnetic 
induction and the related magnetic flux exist only in the 
layer of the thickness A. Meissner currents in this layer 
screen the internal field 4ttM, as well as they screen the 
external magnetic field in a nonmagnetic superconductor. 

Let us consider now the mixed state of the supercon- 
ducting ferromagnetic, in which vortices (magnetic flux- 
ons) are present in the bulk. Since ferromagnetism does 
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not affect the London equation (||) , one expect the same 
magnetic-induction distribution in the mixed state as for 
nonmagnetic type II superconductors and the free 
energy is given by 



F m (M,B) = f E +K- 



2irM 2 



B-M + F (B), (6) 



where Fq(B) is the free energy of a nonmagnetic type 
II superconductor, and B now is the magnetic induction 
averaged over the vortex-array cell. The energy Fq(B) 
contains both the magnetic energy B 2 /8tt and the kinetic 
energy of the superconducting currents inside the vortex 
cell. Determining the magnetic field 



H = 4tt 
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we see that the magnetization curve of a superconduct- 
ing ferromagnet is described by B = Bq(\H + 4ttM\) 
where Bq(H) is the equilibrium magnetization curve for 
a nonmagnetic type II superconductor [7| (Fig. ||a). Note 
that in this relation the magnetic field H has a different 
physical meaning from that used in the Meissner state. 
For the Meissner state we introduced H = B — AirM, 
where the moment M originates from "molecular" cur- 
rents responsible for ferromagnetism, the superconduct- 
ing currents being treated as external currents. In the 
mixed state, which is considered now, it is more conve- 
nient to define the magnetic field as H = B—4w(M+M s ), 
i.e. the definition includes also the diamagnetic moment 
M s = (Bo — H) /Air of the superconducting currents cir- 
culating around vortex lines in the mixed state. Thus 
these currents are treated in the same manner as molec- 
ular currents responsible for ferromagnetism. 

Figure |^b shows that in a superconducting ferromagnet 
the Meissner state (B = 0) exists until H + AirM < H cl , 
where H c \ = ($o/A 2 ) ln(A/£) is the lower critical field 
in a nonmagnetic superconductor and £ is the coher- 
ence length, which determines the vortex core size. So 
ferromagnetism decreases the lower critical field H& = 
H c i — AirM. If AirM > H c \, the Meissner state is absent 
(Fig. ^p) and the superconducting ferromagnet is in the 
mixed state with vortices penetrating into it even in zero 
external field H — 0. This is spontaneous vortex phase 
with nonzero magnetic induction B = Bq(AttM) in the 
bulk. 

Now let us consider formation of the domain structure 
in the standard geometry ||: a slab of the thickness d 
along the anisotropy easy axis y and infinite in direc- 
tions of the axes x and z (Fig. ||). We start from a nor- 
mal ferromagnet. In the absence of an external magnetic 
field the average magnetic induction inside the slab must 
vanish. Therefore, B — in a single-domain structure 
(Fig. ^a), and there exists an uniform magnetostatic field 
H = —AirM in the entire sample, an analog of the elec- 
trostatic field in a charged plane capacitor. This results 
in a high magnetostatic energy H 2 /Sir ~ M 2 . However, 



the domain structure with period I (I -C d) suppresses 
this energy in the domain bulk: H ss and B = —AirM, 
except for the area ~ I 2 near the sample boundary (Fig. 
|b). But the average induction still vanishes, since M 
changes its sign from a domain to a domain. For the 
stripe structure one can solve the equations of magneto- 
statics, Vxff = and V • H = Airp M , exactly [p|,^0|. 
Here pu = —V • M is the magnetic charge. The magne- 
tostatic energy per unit volume of the slab is 

E s = 0.852M 2 ? 2 x i- = 0.852 Af 2 -5- . (8) 
Id d 

This energy is by a factor l/d less than the magneto- 
static energy in a single-domain structure. However, the 
domain walls increase the energy. The energy of one do- 
main wall (per unit length along the slab) is aKSd, where 
6 is the wall thickness and the numerical factor a depends 
on the detailed definition of K and S. Its specification is 
not essential for the present analysis. The domain- wall 
energy per unit volume of the sample is 
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The equilibrium value of the period / is determined by 
minimization of the energy E w + E s 0] : 



I = 



aK 



0.852Af 2 
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Let us return back to a superconducting ferromagnet. 
In the Meissner state the magnetic induction must van- 
ish in the bulk, which is compatible only with the single- 
domain structure. Thus the equilibrium domain struc- 
ture is impossible in the Meissner state. However, do- 
mains with the changing direction of M can appear in the 
spontaneous vortex phase with nonzero B — Bo(AirM). 
Like in a normal ferromagnet, the magnetic flux oc B in 
domains should produce the magnetostatic fields in the 
area ~ I 2 , but in a superconducting ferromagnet these 
fields are by the factor B (AirM )/AttM smaller. We can 
take it into account introducing the effective magneti- 
zation M = B (AttM)/Att. Then the period of the do- 
main structure is given by Eq. (|To|), where M must be 
replaced with M. In the limit of large AirM ^> H c i, 
one has M — > M and the effect of superconductivity on 
the domain structure vanishes. In the opposite limit of 
small M, when 47rA^f — > H c i, M vanishes and the pe- 
riod I becomes infinite, as it should be in the Meissner 
state 47rAf < H c i. However, this calculation of / assumes 
that the penetration of the magnetostatic field into a su- 
perconducting ferromagnet is similar to the penetration 
into a normal ferromagnet. The assumption is correct 
if rigidity of the vortex array is negligible and the effec- 
tive penetration depth is infinite. We can also consider 
the opposite limit of a very rigid vortex array, when the 
magnetostatic fields penetrate only into the layer of the 
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thickness A. If A <C I, the penetration of the magnetic 
flux into a superconductor becomes insignificant. This 
increases the magnetic fields outside the sample, as well 
as the total magnetostatic energy, by a factor of 2, while 
the correspondingly period I decreases by a factor of \/2 
(Fig. ||c) , in analogy with the effect of a superconducting 
substrate on a domain size in a ferromagnetic slab JlO| . 
Thus avoiding a detailed analysis of the vortex and field 
pattern in the domains close to the sample border we lose 
only a numerical factor of not more than In any case, 
superconductivity, which coexists with ferromagnetism in 
the same bulk, always increases the domain size, in con- 
trast to the superconductor-ferromagnet bilayer, where 
superconductivity shrinks ferromagnetic domains fljcfl . 

The absence of the equilibrium domain structure in 
the Meissner state does not rule out a possibility of 
metastable domain walls, as topologically stable planar 
defects. Domains can appear also because of disorder, or 
grain structure. The structure of the domain wall should 
be found by solution of the coupled equations of mag- 
netostatics and the London electrodynamics. We restrict 
ourselves to the simplest case, when the London penetra- 
tion depth A essentially exceeds the domain wall thick- 
ness S. This means that at the spatial scales of order 
6 the domain-wall structure is governed by large ener- 
gies [the exchange energy and the anisotropy energy, see 
Eq. (|j)] and is not affected by the magnetostatic and ki- 
netic energy. On the other hand, at scales ~ A one can 
find the distribution of B and H from the London equa- 
tion at constant M. This is shown for the Bloch domain 
wall (the magnetization M rotates in the plane of the 
wall and does not produce the magnetostatic charges) 
in Fig. [|a. Though our picture corresponds to a 180° 
wall, a similar picture is expected for any domain wall. 
The jump of the tangential component of the moment 
M at the wall defines the current sheet, responsible for 
a jump of the magnetic induction parallel to the wall, 
whereas a possible jump of the normal component of M 
(a "charged" domain wall) would produce a jump of the 
normal component of the field H. 

The magnetic flux on the opposite sides from the do- 
main wall creates the magnetostatic fields outside the 
sample, where the wall meets the sample boundary (Fig. 
^b). The magnetic fluxes, which exit from the sample at 
two sides from the wall, are equal in magnitude (47rMA 
per unit length along the wall) but opposite in direction. 
The magnetostatic field from domain walls could be used 
for their experimental detection. At distances r ^> A 
from a line, where the wall exits to the sample boundary, 
this field is a dipole field of the order of MX 2 /r 2 . 

In summary, the letter presents the first analysis of 
the domain structure in superconducting ferromagnets. 
There is no equilibrium domain structure in a supercon- 
ducting ferromagnet in the Meissner state. In the spon- 
taneous vortex phase the period of the domain structure 



may essentially exceed that in the normal ferromagnet. 
But metastable domain walls can exist even in the Meiss- 
ner state. They generate the magnetic flux in layers of 
a thickness A, which can be revealed by magneto-optical 
methodic. 
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FIG. 1. Magnetic induction B (solid line), magnetic field H 
(dashed line), and AttM (dotted line) at the boundary between 
a superconducting ferromagnet (a; < 0) and vacuum (x > 0). 
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FIG. 2. Magnetization curve: a) nonmagnetic type-II su- 
perconductor; b) superconducting ferromagnet, 4nM < H c \; 
c) superconducting ferromagnet, 4nM > H c \. 
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FIG. 3. Domain structure in normal and superconducting 
ferromagnets. The thick arrows show directions of the mag- 
netic moment M, the thin lines with arrows are force lines of 
the magnetostatic field H. The magnetic charges are shown 
by + and -. a) A single-domain structure. In the whole 
bulk B = and H — —4nM. b) A stripe domain struc- 
ture in a normal ferromagnet. The magnetostatic fields are 
present in areas ~ I 2 inside and outside the sample. In the 
rest parts of domains H = and B — 4nM . c) A supercon- 
ducting ferromagnet in the spontaneous vortex phase with a 
rigid vortex array. The magnetostatic fields appear only in 
areas ~ I 2 outside the sample. In the bulk of domains H = 
and B = B (4ttM). 
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FIG. 4. Domain wall in the Meissner state: a) Magnetic 
induction B (solid line), magnetic field H (dashed line), and 
47rM (dotted line) near the domain wall in the superconduct- 
ing ferromagnet. b) Magnetic flux lines around the exit of the 
domain wall (of thickness 8) to the sample surface. 
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